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Abstract 

In this paper we present a battery of results related to how Galerkin semidis- 



' cretization in space affects some formulations of wave scattering and propagation 

^ ■ problems when retarded boundary integral equations are used. 



1 Introduction 



1> \ In this paper we present several results concerning the effect of Galerkin semidiscretization 
^ ' in space when applied to formulations where time domain boundary integral equations 
OO ■ are used. The focus will be set on energy conservation properties (or lack thereof) and 
^ the model equation will be the acoustic wave equation in two and three dimensions. 

While frequency domain boundary integral equations are very often used for scattering 
problems and for construction of absorbing boundary conditions, their time domain coun- 
terparts have enjoyed a more limited success. The work of Alain Bamberger and Toung 
Ha-Duong [3]-|l] in the mid-80s sparked theoretical and practical activity concerning 
^ ■ time domain integral equations for acoustics, elastodynamics and electromagnetism. A 
detailed account of this early work can be found in the review articles [Hj and [H]. It 
is nowadays well understood that some of the apparently unsur mount able difficulties in 
the development stable discretizations of the time domain integral equations are related 
to the correct approximation of integrals and that these equations can be advantageously 
used in situations where frequency domain formulations are not competitive (scattering 
of frequency rich noise) or applicable (coupling with nonlinear problems). In this paper 
we will study semidiscretization in space with general Galerkin schemes, which can be 
complemented either with Galerkin in time discretization (following the philosophy of 
[3] and its consequences) or with Convolution Quadrature in time (following the seminal 
work of Christian Lubich [2Tj and its development in [5], |H], [24], [7]). 

From the point of view of the type of problems we are going to study, we can classify 
them in two groups: 
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(a) Scattering of acoustic waves by impenetrable obstacles. This leads to an exterior 
wave equation with Dirichlet or Neumann boundary conditions on the surface of the 
obstacles. Instead of dealing with the traditional (and more practical) problem of an 
incident plane or spherical wave (that has no finite energy), we will study the effect 
of the scattering of a free wave that propagates from some compactly supported 
initial conditions. 

(b) Propagation of waves originated in a locally non-homogeneous medium surrounded 
by a homogeneous medium. In this case, we will also follow the propagation of com- 
pactly supported initial data when we surround the computational domain by an 
artificial boundary where we impose non-local boundary conditions using time do- 
main integral equations that give an exact expression of the corresponding Steklov- 
Poincare operators. 

In order to present the techniques and results in a clear way we will first focus on the 
problem of scattering of waves by a sound-soft obstacle, formulated using a single-layer 
potential representation (i.e., with an indirect boundary integral method). We will show 
that any Galerkin semidiscretization in space of the corresponding retarded integral equa- 
tion is stable from the point of view that a certain energy functional is preserved over 
time. In fact, we will show that the total energy that is conserved has to be computed 
counting potential and kinetic energy in the exterior and interior of the obstacle, thus 
giving theoretical evidence of the fact that energy is leaked to the interior of the scat- 
tcrcr, but the total balance of energy is still correct. In a second step, we introduce the 
second class of problems where a locally non-homogeneous medium is surrounded with an 
artificial boundary (that also surrounds the support of the initial data) and we follow the 
propagation of these initial data when we discretize in space with a Galerkin scheme both 
for the variational formulation of the wave equation in the computational domain and for 
a system of two retarded boundary integral equations that are coupled with the interior 
wave equation to construct an absorbing boundary condition on the artificial interface. 
At the discrete level, this can be understood as a time domain coupling of Boundary and 
Finite Elements (BEM-FEM). The result that we show is very similar to the one for the 
scattering problem: energy is conserved as a function of time (independently of Galerkin 
discretization) as long as we count the energy of a ghost wave in the computational do- 
main -the energy of this interior wave has to be computed using the material properties 
of the surrounding medium and is added to the natural energy of the total wave in the 
interior and exterior domains. 

Once these two model situations have been analyzed, we present two very simple ex- 
tensions where the same techniques apply to establish similar results. The first extension 
deals with the scattering of a sound-hard obstacle formulated with a double-layer acoustic 
potential. The second one is the construction of a tighter integral absorbing boundary 
condition using two integral equations but only one boundary unknown. We next show 
that there are formulations where energy is not conserved after space discretization. Two 
of them are explored: the use of a direct integral formulation (based on Kirchhoff's for- 
mula) for the sound-soft scattering problem and the construction of absorbing boundary 
conditions with only one integral equation. While the negative results do not imply that 
these discretizations are unstable, they might indicate some undesirable effects. Actually, 
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the direct integral formulation for the sound-soft scattering problem, fully discretized with 
a Galerkin method in space and Convolution Quadrature (CQ) in time, can be analyzed 
[6] (See below for more on this.) As for the one-equation absorbing boundary condition, 
at present time there seems to be insufficient theoretical and practical evidence to sug- 
gest convergence or lack thereof. The extension of all the results in this paper to linear 
elastic waves is straightforward. Extension to other types of elastic waves (viscoelastic or 
poroelastic, for instance) and to electromagnetic waves will possibly require more effort. 

Apart from the inherent interest of having theoretical evidence of energy preserva- 
tion after space discretization of some integral and integro-differential model problems for 
which this energy conservation property does not seem to be evident, the type of results 
that we prove in this paper has some other applications. On the one hand, the transfor- 
mation of the semidiscrete equations to a transmission problem is allowing us to provide 
an analysis of the CQ-BEM for some scattering problems [6]. Up to the present time, 
analysis of discrete methods for retarded integral equations was based on estimates in the 
Laplace domain. The dynamical system approach that we develop in this paper can be 
used to obtain better estimates of the effect of space discretization and approximation of 
data in some retarded integral equations. It is also our belief that these techniques will 
be advantageous for some kind of analysis of fully discrete BEM-FEM schemes using CQ. 
In this sense, this paper can be taken as the time domain counterpart of what was done 
in [20] in the Laplace domain and we hope to be able to build on these two articles for 
the complete analysis of the fully discrete method. 

This paper offers a novel approach to the analysis of time-domain integral equations 
and their semidiscretization in space. Theoretical study of retarded integral equations 
is almost exclusively carried out in the Laplace domain. Invertibility of the correspond- 
ing boundary integral equations is established in the resolvent set of the Laplacian and 
bounds depending on the Laplace transformed parameter are used to move back to the 
time domain. While this proved to be an extremely powerful approach, leading to solid 
estimates and useful in the analysis of CQ-type time discretization, it is also clear that 
some properties are lost in the process of transforming back and forth to the Laplace 
domain and the some of the functional spaces that this kind of analysis impose might not 
be optimal. This paper offers a complementary point of view using elementary techniques 
from the theory of evolution equations based on strongly continuous semigroups (cf. [12], 
|18j). The main idea is recasting the semidiscrete problem as a Cauchy problem for an 
equation of the second order associated to an unbounded operator. Once the problem is 
reexpressed as a second order equation in time, we use basic estimates in the simplest 
possible setting (Hilbert space theory and a coincidence of two of the four possible spaces 
that appear in the abstract formulation) in order to simply derive the energy conservation 
property. The kind of unbounded operator that will appear in this formulation is very 
closely related to the exotic transmission problems that Antonio Laliena and the author of 
this paper devised for the analysis of CQ schemes applied to non-trivial acoustic scatter- 
ing problems in [2D|. The surprisingly simple (and quite effective) idea consists of dealing 
with the different aspects of Galerkin discretization (the fact that the unknown is sought 
in a discrete space and the fact that the corresponding equation is tested with the same 
space) as non-standard transmission conditions. 

The paper is organized as follows. In Section |2] we present the indirect indirect formula- 
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tion of the sound-soft scattering problem and its Galerkin semidiscretization, and we state 
the energy conservation result for this problem. The proof of this result is carried out in 
Section |3] by using the techniques that have been described before, and taking advantage 
of well known results for the single layer potential associated to the Yukawa (reaction- 
difussion) equation. In Section H] we present a boundary integral absorbing boundary 
condition for a wave propagation problem in free space and its Galerkin (BEM-FEM) 
discretization, and we state the corresponding energy estimate. The proof of this result is 
the content of Section [51 Section |6] presents two easy generalizations: indirect formulation 
for sound-hard scattering and a different absorbing boundary condition. In Section [7] we 
give two partially negative results, showing that space discretization does not lead to en- 
ergy conservation either in a direct formulation of a scattering problem or a one-equation 
absorbing boundary condition. Finally, Section |8] includes the abstract frame of Cauchy 
problems for equations of the second order associated to a class of unbounded operators. 
These results are easy modifications of theorems that are well known in the literature of 
Cq— semigroups and are presented here in order to have an easy reference in the same 
language that we are using in the paper. 



Notational foreword. Basic elementary results on Sobolev spaces will be assumed 
throughout without specific reference. All of them can be found in any text on the 
subject (P for instance). Given an open set U with Lipschitz boundary, we will consider 
the spaces L'^{U) and H^{U) for m > 1. The L^(^)— norm will be denoted || • \\u and 
the H^{U)—Yioim. || ■ \\i^u- Two fractional Sobolev spaces will be used on the boundary 
H^^/\dU). The space H^{U) is the kernel of the trace operator H\U) H^I^{dU). 
The characteristic function of U will be denoted xu- 

Although the geometric layout of the different parts of this article will vary, in all cases, 
there will be a bounded open set fi, with exterior fi"*" := R'' \ and common Lipschitz 
boundary F := dVt = dQ~^. There are two possible traces on F, which will be respectively 
denoted 7+ : H^{n+) H^/'^{T) and 7" : H^{n) H^/'^iT). Whenever there is a single 
trace or functions are exclusively defined in the interior domain, the superscript will be 
dropped. 

For functions of space and time variables m(x, t), we will often employ the notation of 
theory of evolution equations, where only the time variable is displayed. This amounts to 
considering functions u : [0,T] — > X, where X is a space of functions of the x variable. 
With this notation, u and il denote the first and second t— derivatives of u{t). 



2 Single layer potentials for sound-soft scattering 

Consider a bounded open set f2 C M'', with Lipschitz boundary F := dQ such that 
Q'^ := M'^ \ is connected (f2 does not need to be connected though). Let uq and Vq be 
given functions such that O := supp Uq U supp Vq is compact and O HQ = 0. We consider 
the following problem of scattering of acoustic waves by a sound-soft obstacle: we look 
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for u : [0, oo) H^{n+) such that for alH > 

u = Au in , (la) 

u = on r, (lb) 

and the initial conditions (on Q~^) 

u{0) = uo, u{0) = vo (2) 

are satisfied. We will additionally assume that Uq G HQ{fl~^) and Vq G Hq{Q^). Let R > 
be such that 

Hue C 5(0;i?) := {x G R'^ : |x| < i?}. (3) 

Because of the finite speed of propagation of waves, solutions of problem (II])- (12]) are 
supported in 5(0; R + T) for all t G [0, T]. 

The solution of ([I])-([2]) will be decomposed as the sum of a free (or incident) wave and 
a scattered wave. To define the free wave we need to extend the initial data by zero to 
the interior of the scatterer: 

( n " o'' 

m \ l. 

The free wave is the solution of the wave equation in R'^ for t > 

ii^'^^ = Au^'^^ (4) 

satisfying initial conditions 

„fr-(0) = Euo, ii^'^'iO) = Evo. (5) 
The solution to this problem satisfies 



u 



free 



G C2([0, oo); L2(R'^)) n C\[0, oo); H\R'')) n C([0, oo), //'(R'^)) (6) 
and supp tt^'''"^(t) C B{0] R + t) for all t. The scattered wave is the difference u — u^^'^'^. 

Remark 2.1. For the kind of arguments that we are going to develop, we will assume 
that the free wave is known. Section HI deals with a semidiscrete version of a problem 
that generalizes dl])-©, leading to algorithms to compute the propagation of compactly 
supported initial data in free space. In applications, the free or incident wave is actually 
known. Typically the free wave is the solution of a non-homogeneous wave equation and 
has singularities at some points away from the scatterer. Other practical incident waves 
include plane waves that are not compactly supported at any given time and have infinite 
energy. Because our interest lies in how discretization affects the integral model used to 
represent u — u^'^'^'^, the use of free waves as those of (ll])-([5]) fulfills our needs. 

The scattered wave u — u^^'^^ can be represented with a single layer retarded potential 
[3J. For (x, t) G f2+ x [0, oo) and a given density i/j : [0, oo) — )■ H~^/'^(T) the single layer 
potential is defined as 

f V^(y t -|x-y|) (whenrf = 3), 

(5*^)(x,t) = <( . /.t-|x-y| ,/;rv 

/ / ' dr(y)Hr (whenrf = 2). 

hk V(t-r)2-|x-y|2 ^ ^ 
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The integral forms in ([7]) are only valid for smooth densities. Weak forms of the potentials 
have to be used in general (see [19] for full justification in the three dimensional case). We 
will keep the convolutional notation for the layer potential 5 * -0 in order to distinguish 
time domain potentials and operators from similar entities for steady state problems. Let 

V * ip := 7(5 * ip) 

be the corresponding single layer integral operator (see [3] and ^^)■ The indirect repre- 
sentation of the scattered field looks for a density ip -.M ^ H~^^'^{T) such that 

= in (-00,0) V * ^ + ju^""" = (8) 

(this is an equation on F x [0, 00)) and then represents the total wave by 

u = S*ij + u^'^\ (9) 

Equation ([8]) is now approximated with a Galerkin method only in the space variable. To 
do that, we choose a sequence of finite dimensional spaces 

approximate ([8]) by the problem 

iph-.^^ Xh, ^/'/i = in (-00, 0), 

(10) 

and write an approximation of the total wave as 

Un:=S* + (11) 

The angled bracket in ([10]) denotes the H-'^/'^{T) x H^/'^{T) duality product. Note that 
suppM/i(t) C -8(0; R + t) for all t, because the layer potential S * ip propagates a wave 
from r at unit velocity starting at time t = 0. 

Remark 2.2. Problem (1101) is a system of functional equations. Let us assume that 
Xh C L°°{T) and that {Nj : j = l,...m} is a basis of X^. We can then write the 
unknown density as 

M 

My, t) = Yl W : K ^ M, supp C (0, 00). 

i=i 

In the three dimensional case, ([TOl) is then equivalent to 

E / / |x-y|)dr(x)dr(y) = - f iv,(x)txf-^(x,t)dr(x) 

i = l,...,M. (12) 
This system of functional equations include integrated delays of all the unknowns. 
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We are now in conditions to state the main theorem of this part of the article, showing 
energy conservation and regularity in time for the semidiscrete total wave field. The proof 
of this theorem will be given in Section |3l 

Theorem 2.1. The semidiscrete total field Uh and the associated density i/jh, given by 
(^-(^, satisfy 

Uh e C\[0,oo);L^iR'^'))nC\[0,oo);H\R'^)), (13) 
G C{[0,ooy,H-'/\T)). (14) 

Moreover, the energy 

is constant over time. 

Let us emphasize that the exact total wave field u is only defined in and that it can 
be extended by zero to the interior of the obstacle. The semidiscrete total wave field Uh is 
however defined in M'^ and the property of energy conservation is proved with integration 
over and not only over fl~^. This shows that Galerkin discretization leaks part of the 
energy to the interior of the obstacle, although the total energy is still constant. 



3 Proof of Theorem 12.1 

Notation. The jump of the trace across F is denoted 

[7?/] := — 7^M. 

Weak normal derivatives on F can be defined using Green's formula. Given an open 
bounded set B that contains Q, we can define d^u for any u G /f^(B \ F) such that 
Au G L2(B \ F) with the formulas: 

{d-u,j-v)r = iVu,Vv)n- + iAu,v)n~ e H\n-), 

(a>,7+z;)r = -iVu,Vv)Mnn+-iAu,v)Mnn+ Wv e H^^iM D n-^) , 

where 

H^^{M n ^]+) = {veH\Mn n+) ^qmv = 0} (I6) 

and 7aB is the trace operator associated to the boundary of B. It is well known that the 
definition of the exterior normal derivative is independent of the set B. The jump of the 
normal derivative [d^u] := S^-u — d'^u is defined for the same class of functions. 

The single layer Yukawa potential. Consider the fundamental solution of the Yukawa 
operator u i— ^ — Au + u: 

■Ix-yl 

(when d = 3), 



E(x,y):=<^ 47r|x-y| (17) 
-Ko(|x-y|) (whenrf = 2), 



27r 
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where Kq is the modified Bessel function of the second kind and order zero. On the 
surface/curve E := dMU T, we can define the single layer potential 



SA :=^i?(.,y)A(y)dH(y). 



Using a weak definition of this potential (see [TU] or the very general theory developed in 
[22]), we can prove that S : H-^/^{E) H^R'^) is bounded, that 

V := 7hS : R-^/^E) H^/\E) (19) 

is bounded and coercive 

(A, VA)h > VA G H-'I\E), (20) 

and that A(SA) = SA in M*^ \ S for all A. Also, if we write A = (Ar, Xa) e H-^/^{E) = 
i7-i/2(r) X H-^/\dM), then |a^(SA)] = Ar (the jump is defined across T). 

A transmission problem. The single layer retarded potential satisfies the following 
properties [3] 

[7(5*^)1 = [9.(5*^)1 = ^ (5*V')(0) = |(5*V')(0) = 0. (21) 

For given initial data with compact support we can choose i? > so that ([3]) is satisfied. 
Since the speed of propagation of u^'^'^^ and S * i/jh is the same, 

suppuhit) CM := B{0;R + T) \/te[0,T]. (22) 

Noticing that 

[7« J = tiS * ^P,)] + = 0, (23) 

it follows that, restricted to the time interval [0,T], the function defined by (fTO|) - (fTT|) can 
be understood as Uh [0,T] ^ H^{E \ T) that solves the wave propagation problem: 

Uh = ^±Uh, (24a) 

buh\ = 0, (24b) 

e XI (24c) 

[d^u^] e A,, (24d) 

7dBUh = 0, (24e) 

with initial conditions 

Uh{Q) = Euo Mh(0) = Evo. (25) 
The set in fl24cp is the polar set of Xh, i.e., 

A° := {e e H"\V) : (/i^,e)r = V/i;, G Xu}. 

The Laplace operator A± in fl24ap is the one defined in the sense of distributions in B \ F. 
Finally, the transmission condition (I24d|) is equivalent to 

([a.«J,Ur = o 'iihexi 

Let now Uh be a solution of (!24|) - (l25|) and define ifjh '■= \d,jUh\ = [duiuh — u^^'^'^)} (recall 
) . We can then show that S * iph = Uh — u^^'^^ by comparing the transmission problems 
that both solutions satisfy. 
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Formulation as a Cauchy problem. Consider the Hilbert spaces 

H := L^(E), 

V := {ue ifo'W : 7U e X°}, 
D{A) := {ueV : A±u G ^^(B), {d^uj G X^}, 

endowed with the respective norms 

1 /2 

||«||// := Iklll := ||Vu||b, := (||Vm||b + ||A±m||^) 

We also consider the operator A := A±. We next verify the two conditions of Section [HI 
First of all, the generalized Green's Identity: for u G D{A), v E V C Hq(E), using the 
weak definition of the normal derivatives, it follows that 

(Vii, Vv)m + (A±M, ^;)b = {ld^u},'yv)r = 0, 

because [duu] G Xh and 'yv G X'f^. The second step is surjectivity: for any / G L^(B) we 
want to find 

ueD{A), -Au + u = f inl\r. (26) 

We first choose G H^{M) such that -Am^^ + = / in B, and note that m"^ G H^{M) 
by a simple regularity argument. We next consider a variational problem in the space 
X^ := Xh X H-'/\dM) C //-i/2(H): 

A = {Xh, Xd) G Xh^ 

27 

. (P, VA)h = -(p.,7«°')r Vp = (p,,pa) G X,. 

This problem is uniquely solvable by the coercivity property (!20|) . We finally take u := 
u^'^ + SX. It is clear that u G i/^B). Testing ([27D with elements (0, po) G {0} x //-^/^(SB) 
and recalling that 7=8 = V it follows that u G Hq{M). Testing with elements (p/i,0) G 
X/, X {0} it follows that G X'^. This proves that u eV. Also, -A-u + -u = /inB\r 
and [duu] = [duSfi] = A/i G X/^. Therefore u G -D(A) and (|26l) is satisfied. 

Conclusions. The theory for abstract Cauchy problems that is sketched in Section 
[H] can be applied to problem (^^. When initial data are in D{A) x V, the Cauchy 
problem has a unique strong solution. Since Uq G H'^{fl~^) and G vanish 
in a neighborhood of Q and also in a neighborhood of 9B, then (Euo^Evq) satisfies the 
regularity requirements to have strong solutions: 

Uh G C\[0, T]; L\M)) n C\[0, T]; H',(M)) H C([0, T]; 

Extending by zero to the exterior of B, it is clear that f fTSj) is satisfied. Note that the 
property is proved in [0, T] for any T and that T influences the size of the ball B that we 
use to truncate the domain. Note also that ■ ] : D{A) — )■ H~^^'^(T) is bounded. Then 

A-= Id.Uh] eCi[0, 00); H-y^T)). 

Remark 3.1. The convenience of adding the cut-off boundary dM -far enough from the 
obstacle so that the solution is not affected in the time interval [0,T]- is due to the 
difficulty of working with energy spaces in unbounded domains (see Section [8]). 
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4 Transparent boundary conditions 

We now consider the problem of propagation of compactly supported initial conditions in 
free space M*^ 

c'^u = V ■ (kVu) 

with initial conditions 

U{0) = Uq, u{0) = Vq. 

We assume that k,c & L°°(]R'^) are positive and that k^^,c^^ G L°°{ 
we assume that c = 1 and k = 1 outside a compact set. Let 



). Furthermore, 



O := SUpp Uq U SUpp Vq U SUpp (c — 1) U SUpp {k — 1) , 

which is a compact set by all the previous hypotheses on initial data and coefficients. We 
choose a bounded open set Q with Lipschitz boundary F := dQ, with connected exterior 
f2+ := M'^ \ F and such that O C Vl. We admit the possibility of taking n = O, as long 
as O meets the regularity hypotheses required for Q. Also, there is no need to have Q 
connected, although for simplicity we will assume that its exterior is connected. 



Denoting the conormal derivative on F by d^u = (kVu) ■ u, and renaming u' 



ext 



we have the problem 



u 



ext 



in Vl , 
in n+, 
on F, 
on F„ 



with initial conditions 

m(0) = Uq, {i(0) = Vq, 1x'=^*(0) = 0, M°^'(0) = 0. 

The function u^^^ can be represented with Kirchhoff 's formula 



,ext 



U 



V*ip-S*X iip,X):= (7+«^^\ d+u^'"'), 



u\n+, 

(28a) 
(28b) 
(28c) 
(28d) 

(29) 

(30) 



using both Cauchy data on the interface F. This formula employs the double layer retarded 
potential (see [i]. p^) 



y|) 



r V 47r|x-y| 



z=y 



^(y)dF(y) 



(when d = 3) 




r Jo 



^ _ _ ^y^i^,y,t - T)dT{y)dT 



y|). 



where 



'f(x,y,t) 



-$(x,y,t-r)dF(y) 
(x - y) ■ jy(y) 



(when d = 2), 
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and ^'(y) is the outwards pointing normal vector at y. Three retarded integral operators 
appear in this formulation: 

)C'*X:= \{d+ + d;){S*X) /C*^= i(7+ + 7-)(I)*(^) W*^ = -d^{V*^). {31) 

The boundary- field formulation is based on representing u^^^ with the formula fl30|) . We 
then write a weak formulation for the interior equation ( ]28aP and substitute the trans- 
mission condition (!28dll : 

{c-%v)n + {KVu,Vv)n- {X,iv)r = Wv e H\fl). (32) 

A second equation is obtained by imposing the transmission condition f l28cp . with u^^^ 
written in terms of its Cauchy data fl30l) and using the conditions (|3T1) to represent the 
trace of "D * 

7u + V*A-(iy? + /C*^) = 0. (33) 
The third equation is an identity satisfied by the Cauchy data: 

iA + /C* * A + W * = 0. (34) 

Finally, the global formulation results from writing ( l32l) -(l33l)-(1Ml). The unknowns are 

with initial conditions 

^(0) = uq, u{0) = Vq, a = in (— oo, 0), 99 = in (—00, 0). 
We now choose spaces 

Vh C H\n) Xh C H-^'\V) Yh C H^'\V), 

and look for 

Uh : [0, 00) y/,, Xh-.R^ Xh, V^h : M ^ Yh, 
satisfying initial conditions 

M/j(0) = Mfe,o, Uh{^) = Vhfl-: A/i = in (-00,0), v'/j = in (-00, 0), (35) 

for approximations uq ^ Uhfi G Vh and vq ~ Vhfi G Vh to be determined by a projection 
method or by some kind of interpolation process (this is not relevant in the sequel). 
Finally, we have the set of Galerkin equations: for t > 

(c"^iift, Vh)n + {i^'^Uh, Vvh)n - {Xh, lVh)v = Vv,, G Vh-, (36a) 
(;Uh,7Mh)r + (/U,,, V * A/,)r - (/U/i, Iv^h + /C * V9h)r = Wfih & Xh, (36b) 
qXh + K} * A,,, ih)v + (W * Vh, a)r = Va G (36c) 

Equations (136!) form a system of linear second order differential equations coupled with 
the kind of retarded equations that we found in Section [3] (see (fT2l) for instance). These 
equations are complemented with the initial conditions ( l35l) . Once (!35l) - (l36l) has been 
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solved (its solvability is part of what we state in the next theorem), we can define the 
approximation to the exterior solution 

u'^' := V * iph - S * Xh. (37) 

Note that this function is defined in M*^ and not only in Q~^. In order to simplify some of 
the forthcoming arguments, we will assume that constant functions belong to the three 
discrete spaces 

Po(^^) c Vh Po(r) c Po(r) c y,,. 

When r and Q are not connected, the spaces fo(T) and Po(f^) have to be understood 
as the spaces of constant functions on each connected component of the corresponding 
domain. 

The following result (which will be proved in Section E]) gives a basic regularity estimate 
for the solution of this problem. It also states an energy conservation property, where in 
addition to the expected wave fields {uh in the interior domain Q and m^^* in the exterior 
domain we have to count the energy of u'j^^ in the interior of F, computed with the 
material properties of the surrounding medium. 

Theorem 4.1. The semidiscrete total fields {uh,u'f^^) and the approximations to the 
Cauchy data on the interface {iph,Xh) given by fl35|) - fl36|) - fl371) satisfy 



Uh e CWO,oo);L\n))r]C\[0,oo);H\n)), (38) 

Mf* G C\[0,oo);L^{R'^)nC\[0,oo);H\R'^\T)), (39) 

A;, G C{[0,ooy,H-'/\r)), (40) 

iPh G Ci([0,oo);ifV2(r)). (41) 

Moreover, the energy 



is constant over time. 



5 Proof of Theorem 4.1 



Introduction of a cut-off boundary. We will prove the theorem for an arbitrary 
interval [0,r]. Since u^"^* is defined with retarded potentials whose densities are causal 
functions (see the initial conditions f l5^ ). we can pick a sufficiently large radius > so 
that 

supp<"*(t) C B := 5(0; R) Wt G [0, T]. (43) 

A transmission problem. The first step towards the proof consists of writing (136|) in 
terms of the fields {uh,u'f^^). The field u'f^^ satisfies the wave equation in M*^ \ T. The 
interior field Uh does not satisfy a differential equation though. Note that: 

[7«ri = -^H [a.^ri = -A. (44) 



12 



The transmission problem looks for 

Uh : [0, oo) Vh, 
satisfying initial conditions 

UhiO) = Uhfl, Uh{0) = Vh,o 
and the equations for alH > 



[0,00) ^ H\M\T) 



(45) 



(46) 



Ur1,7^/^)r = 
0, 



(47a) 
(47b) 
(47c) 
(47d) 
(47e) 

Equation fl47ap corresponds to fl36al) after substituting A^, = — [c^i/tt^''*] (see (jH]). The 
exterior boundary condition fl47bp is a consequence of fj43|) . The first condition in both 
fl47dp and fl47ep is a consequence of (jH]). The second condition in fl47dp is just fl36bl) . 
Finally, the second condition in fl47ep is fl36cp . 

Formulation as a Cauchy problem. Consider the space 

Hl^{M \ T) := {u e H\M \ F) : janu = 0} = H\n) x H'g^iM f] fi+) 
(recall 0160 ). The three required spaces to fit in the frame of Section [S] are: 



H 
V 

D{A) 



(4J 



{{uh, u'') eVhX H'q^{M \ F) : M e Yh ju^ - T+w'^ e Xn (49) 

(50) 



\\iuh,u*)\\j, := \\c'W\\l+\\u*\\l. 



The norm of H is 
The following norm 



{'ju — 7^m) 



+ 



1 

W\ 



u 



can be easily shown to be equivalent to the usual Sobolev norm in this space. This allows 
us to write a norm in V: 



1 



u 



since we have assumed that Po(F) C Xh- For D{A) we define the norm 

IIK,«*)IId(A) := IIK,«'')llv + l|A±«1|^\r- 
13 



To define the operator A : D{A) — )■ H associated to the evolution problem fH7|) . we 
introduce the operators : H^(^l) — )■ Vh and 7^ : H^^^'^{T) — )■ T4, defined by the 
discrete equations 

(c-^Alu, Vh)n = -(/tVw, Vvh)n \fvh E 14 (51) 

and 

{c-^liX,Vh)n = {X,7Vh)r E Vh, (52) 

respectively. These operators are defined so that fl47ap can be rewritten as iih = ^l^h — 
^lXduU^jf'\. The operator A : D{A) ^ H is then defined by 

A^u^y) := {/^W - li{dXlA±u''). (53) 
Problem fH7|) has the general form fl87|) with initial data {uhfl, 0) G D{A) and {vh^, 0) G V . 



Rigid motions of the system. Let M := {0} x spanjxn}- Noting that 

IduXn] = 0, I'jxnj = 1 eYh, 'yanXn = 0, 5^Xn = 0, and = 0, 

it is simple to check that M C D{A) and M C ker(j4). We consider the following 
seminorm in V 

\iuh,u*)\^y := \\K'/^Vuh\g + \\\/u*\\l\^, 

associated to a semi-inner product [{uh,u*), {vh,v*)]v- The hypotheses to consider the 
finite dimensional space M as a space of rigid motions of the evolution problem (IHTI) (see 
Section [HI) are then easily verified. 



Verification of the associated Green's Identity. Let u := {uh,u*) G D{A) and = 
{vh,v*) G V. Using the definition of the discrete operators fl5Tl) - fl52l) and the definition of 
the weak normal derivatives, it follows that 

iAu,v)H + [u,v]v = {c-^Aluh,Vh)n - {c-^-rild^u*},Vh)n + {A±u\v'')M\r 
+ (nWuh, \/vh)n + (VM^ Vt;*)B\r 

= -(|yi, j^;,-7+< )r - {d^,^rifY)T = 0. 



The Yukawa double layer potential. For the proof of the surjectivity we need to 
introduce the double layer potential for the Yukawa operator on S = F U SB 

Bip:= J^VyEi ■ ,y)y.(y)dH(y), 

where E is given in ( fTTj) . This potential defines a bounded operator 

B : H^/^{E) ^ H\R'^\E) (54) 
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such that A(D(/j) = Dtp in M'^yS for all ip. Two bounded integral operators are associated 
to this potential 

K := i(7± + 7^)D : H^/\E) ^ H^^^E), W := -^J^D : H^/\E) ^ R-^^^E). 

The operator W is coercive [22] 

(Wv^, > Vv. G /fi/^(S). (55) 

Finally the adjoint of K satisfies K* = + djJ^)S, where S is the single layer potential 

defined in ([18]). 

Verification of the surjectivity property. Let {fh, /) G x L^(B) = H. As we did 
in Section [31 we start by finding u""^ G H^{M) n H^{M) such that -Aix°^ + = j b. 
We then look for 

UH eVn, A = (A^, Aa) G X;, := x //-^^(^i) c //-^/^(S) 
^={^h,0)eY,:=Y,x{0}cH'/'iE) 

satisfying the equations: 

(c"\/,, Vh)n + {i^^Uh, Vvh)n - {K, 1Vh)r = (c^'^fh, Vh)n Vv/, G Vh, (56a) 
(/Xft,7Mft)r + (/U, VA)5 - |(/? + Kv9)h = (/i/i,7w'''')r V/x = (/i^^, ^g) G X^, (56b) 
(iA + K*A, Ob + {WV, Oe = {d.u-'', Ch)r = (^z,, 0) G F,. (56c) 

Note that because of the particular form of the space Y_i^, the three bilinear forms where 
either ^ or 99 appear are actually duality products in F. The bilinear form of problem f l56|) 
is coercive in H^{^1) x i/~^/^(S) x if^/^(H) by ([20]) and ([55]) . Therefore problem ([56]) has 
a unique solution. The final step is the verification that the pair 

{Uh,u') := {uH,u^'' + B^-SX) 

belongs to D{A) and that (/ — A){uh,u*) = {fh,f)- This follows from several simple 
arguments that we next list. Because of the potential form for u* and the smoothness of 
u"^^ across F it follows that 

ld,u*j = -A, G Xj,, |7iii = G y,. (57) 

Also 

A±u* -u* = f (58) 

and therefore A^u* G Iv^(B), while u* G i7^(B\F) because u^^ G iJ^(B) and the mapping 
properties of potentials ([T^ . ([M[) hold. Substituting ([FH) in (I56ap and using the definitions 
of the discrete operators ([5Tl) - ([52[) it follows that 

uj,-A^U!, + ^ild,u*j = U. (59) 
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If we test ( l56b|) with elements (0,^a) G {0} x H^'^{dn) it follows that 7^^^* = 0. If we 
test with {nh, 0) G Xh x {0}, it follows that 

^Uh - j^u* = lUh - VA - iliph + Kip) + 7^/"*^ E XI 

Finally, equation fl56cp is equivalent to asserting that 

dtv* = d.v^'' - \\h + K*A -W^e Y°. 

The preceding arguments have shown that {uh,u*) E D{A) while (l58l)-(!59l) proves that 
(J-A)K,w*) = (/,,/). 

Conclusion. We can now apply the theory for Cauchy problems exposed in Section [HI 
The statements of Theorem 14.11 are a direct consequence of these results. 

6 Two simple extensions 

We now show two other situations where the techniques developed in the preceding sec- 
tions can be applied with minor modifications. 

6.1 Double layer potentials for sound-hard scattering 

Let us consider again the geometrical setting of Section [21 The problem of sound-hard 
scattering by an obstacle occupying the region can be expressed with the equations 
([T])-(l2]) with the boundary condition flTbl) substituted by d^^u = on F for all t. The 
assumptions on the initial data are the same as those given in Section [2l If we define u^^*^*^ 
with (j4])-([5]), the solution to this problem can be expressed as u = V * Lp + u^'^'^'^, where 
(/? : M — )■ H^^^(T) vanishes identically for negative values of t. 

The semidiscrete formulation follows from choosing a discrete space 

Po(F)cy, C//^/2(F), 

discretizing the boundary condition 

iPh ■■ '^^ Yh, v^/j = in (-cx),0), 

(60) 

. (-W * iph + d,u'-^-\ a)r = Va G Yh Vt, 

and proposing 

Uh = V*ipn + u'-'''' (61) 

as approximation of u. Note that the following transmission conditions are satisfied for 
all t > 0: 

[7?xJ E [a.iij = 0, <9,n, E Y°. 
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Theorem 6.1. The semidiscrete total field Uh and the associated density given by ( 160|) 
and (16T|) satisfy 



Uh 



Moreover, the energy 



is constant over time. 



e C\[0,ooy,L\R'^))nC\[0,oo);H\ 
G C\[0,ooy,H^/\T)). 



\r)), 



(62) 
(63) 

(64) 



Proof. The techniques for the proof of this resuh are very similar to those of Section HI 
We wiU simply sketch the main steps. First of all we pick a fixed time interval [0, T] and 
choose a ball B that contains the support of the solution for all t G [0, T]. It is now more 
convenient to include a Neumann boundary condition on the cut-off boundary dy^g^Uh = 
for all t. The spaces for the formulation as a Cauchy problem are the following: 



H 
V 

D{A) 



{u G H\M\T) : 
{ueV : A±u G 



en}, 



0, d^ue Y^, dy^Qnu 



(65) 
(66) 
0}. (67) 



The relevant norms and seminorms are: 



\u\v 



\Vu\ 



\u\ 



V 



\u 



V + 



1 f 


2 


1 f 


-. — r / u 


+ 











The associated operator is the same as in Section HJ namely A := A±. The space of 
associated rigid motions is two-dimensional M := span{xf2, xenf^+j- The proof of the 
corresponding Green's Identity is straightforward. For surjectivity we proceed in two 
steps. Given / G L^(B), we first choose G H^(E) such that -Au^'^ + u''^ = / in B and 
diy^dMU^^ = 0. (Note that this is possible by basic regularity theorems of elliptic problems 
on smooth domains.) Next we look for (p = {iph,^d) ^ Uh '■= Yh x H^/'^{dM) C if^/^(H) 
that solves the coercive variational problem 



ve = (a,6)ei:„ 



and define u 



u 



nh 



+ D(/?. (Notations for Yukawa potentials and operators are those of 



Section [H) It is simple to prove that u G D{A) and u 
proof. 



Au = f, which finishes the 

□ 



6.2 A tighter transparent boundary condition 

We now consider the transmission problem (!28|) -( |29|) of Section HI Instead of the trans- 
parent boundary condition that uses approximations of both Cauchy data on F it is 
possible to construct another boundary-field formulation in the spirit of the BEM-FEM 
coupling schemes of Martin Costabel |9j and Houde Han [15]. The boundary unknown 
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is A := dpU^^^. The identities of Cauchy data on the boundary and the transmission 
conditions fl28cp -( l28d|) can be used to produce the following equations: 

|7n - /C * 7n + V * A = - a> = W * - iA + /C* * A. (68) 

In (jnSD and all similar future expressions, it will be understood that ^u{t) = for t < 0, 
although u{t) will be defined for t > and we will not look for a smooth continuation of 
u for negative values of t. 

The boundary-field formulation looks for u : [0, oo) — )■ H^(n) and A : M — )■ H^^^'^{T) 
satisfying initial conditions 

m(0) = -uo, 'u(O) = vq, a = in (— oo,0), 

and the equations 



r 



-(|A-/C**A,7t;)r = Wv e H\n), (69a) 
(/i,i7w-/C*7M)r + (/i,V*A)r = V/x G //-^/^(r), (69b) 
for all t. The discrete version of these equations uses two spaces 

Po(fi) cVhC H\n), Po(r) cXhc h-'/\t), 

and looks for : [0, oo) — )• Vh and A/^ : M — )■ such that 

UhiO) = Uhfl, UhiO) = Vh,o, A/, = in (-cx),0), (70) 

and for all t: 

-{lXh-IC'*Xj,,^Vh)r = yvheVh, (71a) 
{lJ-h,llUh-}C*'yUh)r+ {fJ-h^* Xh)r = M^ih^Xh. (71b) 

Compared with the semidiscrete system ( 136|) . this system has one less group of equations 
and unknowns. The price to pay is the fact that Uh is affected by integral delay operators. 
The reconstructed exterior solution is given by 

■=V*^Uh-S*Xh. (72) 

As in Section m we have to consider u'^jf' defined in R'' (and not only in Vt^) to have the 
correct balance of energy. 

Theorem 6.2. The semidiscrete total fields {uh, u'j^^) and the approximation to the normal 
derivative on the interface Xh given by (!70|) - (!7T]) - fl72|) satisfy 

Uh e C\[0,oo);L\Q))nC\[0,oo);H\Q)), (73) 
e C''{[0,oo);L\R'^)nC\[0,oo);H\R'^\T)), (74) 
A, G C{[0,oo);H-'/\r)). (75) 

Moreover, the energy 

l\\K'/'Vuhml + l\\VuT\t)\\l,^, + l\\c-Wml + U^T\t)U^\r (76) 

is constant over time. 
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Proof. The proof has a similar structure as that of Theorem I4.1[ We will just point out 
the milestones of the proof. Choosing a ball B that contains the support of the solution 
for all if: G [0,T] and adding a Dirichlet boundary condition on dM, we can consider a 
Cauchy problem satisfied by the pair {uh,u'j^^) with the following elements: the spaces 



H 
V 

D{A) 
the norms 



{{uh,u*) e H : lju*j+-fUh 
{{Uh,u*) G V : A±M^ G 



0, -y u* e X^, 7aBM* 



0}, 



\\[Uh,u Jii^ 
and the operator 



|c W\\l + \\u*\\l 



[Uh,U 



\K 



(see (EH) and (l52l)). The associated Green's Identity is straightforward to prove. To 
show surjectivity oi I — A, we proceed as follows. Given {fh, f) G H, we first take 



u 



nh 



G i/^ 



satisfying — A-u + w = / in B. We then construct 



where 

Uh e Vh, 

is the solution of 



{uh, u*) = {uh, u""" + B-fUh - SA) 
A = (A,, Xa) G := X H-'/\dM) c //-^/^^S) 



(77) 



(c '^Uh,Vh)n + {KVuh,Vvh)n 
+ {W^h,'mh - (|A - K'X,^hh = {c-^fh,Vhh + (a.M^^7^'/^)^ G F,,, (78a) 

lr>^h - Ktu/,)h + (/U, VA)h = {^ih, 7u''^)r V/i = {fit, fia) e X^, (78b) 

and jUh = {^Uh,0) G H'^/^{T) x {0} C H^''^{~). Note that problem ([7S]) is associated to 
a coercive bilinear form and is therefore uniquely solvable. It is then simple to prove that 
{uh,u*) G D{A) and A{uh,u*) = {fh,f), which finishes the proof. □ 



7 Two negative results 

This last section shows two results where Galerkin discretizations lead to problems with 
non-constant energy. This fact does not mean that the discretizations are not valid (they 
can still be stable), but at least shows how delicate the energy balance is when discretized 
integral operators are used. 



7.1 A direct method for sound-soft scattering 

Let us consider again the problem of Section [2l The solution of ([I])- ([2]) can be decomposed 
using the free wave (solution to (j4l)-(l5])) and Kirchhoff's formula for the scattered wave: 

u = u^''"' - V * 7^^'-'=" - 5 * A, A := <9+M - d^u^''', (79) 
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where A : M — )■ H^^^'^{T) is causal and in the potential expression T) * ju^^^"^ we have to 
understand that ju^^*^^ : M — )■ H^^'^{T) is a causal function even if w^''^'^ is only defined for 
positive values of t. The indirect decomposition ([9]) used u = S * i/j + u^'^'^'^ was naturally 
extended to M*^, with the result that u{t) = in f2 for all t. The extension to the interior 
domain of (17^ is 

u = u^"''xn+ - 'D * -fu^''^'' - 5 * A, 

because the potential expression P * 'yu^^'^'^ + S * X vanishes identically in Q by Kirchhoff 's 
formula. The discrete version of this process computes 

'Xh-.R^Xh A;, = in (-00,0), 

(80) 

. (/ih, V * A^)r = {fih, h^'"'' - * 7^'"'')r ^^ih e X^, Vt, 
and then constructs the total field 

Uh := - V * ju''''' -S*Xh. (81) 

Note that Uh satisfies the non-homogeneous transmission problem 

Uh = A±Uh, (82a) 

buh] = 0, (82b) 

G X°, (82c) 

[d^Uhj + d^u'^'"' e Xh. (82d) 

(Compare with (12^ and note that the condition at dM can always be added for finite 
time intervals.) The fact that this problem is a non-homogeneous version of a problem 
that is conservative gives a first hint that the natural energy of this problem will not be 
constant. Also, applying integration by parts and fl82|) . we can prove that 




This shows that energy is not constant. 



7.2 Transparent conditions with one equation 

Consider again the propagation problem in free space ( l28l) . Instead of using two integral 
identities as in Sections H] we can work with a single integral equation in the spirit of the 
one-equation coupling of BEM-FEM of Claus Johnson and Jean-Claude Nedelec [171 123] . 
After space Galerkin discretization, the coupled system becomes an evolution problem 
that looks for Uh : [0, oo) — )■ Vh and A/i : R — )■ Xh such that 

Uh{0) = Uhfl, UhiO) = Vh,o, A,, = in (-oo, 0), 

and for alH > 

{c~'^Uh, Vh)n + (kVm/j, Vvh) - {Xh, 'jVh)r = Vv^ G Vh, (83a) 
{lJ.h,l7Uh-IC*'yuh)r + {lJ'h,V*Xh)r = ^fih^Xh. (83b) 
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Note that the first of these equations coincides with (I36ap (the first discrete equation 
of the method in Section H]), while the second one is (I71b|) (the second equation in the 
method of Section [6^ . The discrete exterior solution is defined with (1721) . The same kind 
of manipulations that we have been applied above shows that the pair {uh, "U^^*) satisfies 

with the discrete operators defined in fl5T]) - fl52|) . It is simple to see that the transmission 
conditions above lead to a choice of spaces like flB^ -f E7|) . while the operator itself is 
f lS^ . This mismatch between domain of the operator and operator leads to lack of energy 
conservation, namely, for smooth solutions 

Furthermore, by recasting the evolution problem as a first order system, it is possible 
to show that the corresponding operator A in fl89l) is not maximal dissipative (cf. [HI 
Chapter 4]) and therefore it cannot be the infinitesimal generator of a contractive strongly 
continuous semigroup. As in the case of the discrete Kirchhoff formula of Section ITTTl this 
does not mean that the discretization leads to an unstable method, but it is at least a 
hint that some problems might arise due to the lack of energy conservation. 



8 Some abstract arguments about wave equations 

In this section we summarize the abstract results on Cauchy problems for second order 
equations that we have used throughout the article. The results are elementary conse- 
quences of the theory of strongly continuous groups of isometrics in the Hilbert setting 
that can be found in basic texts as [12j or [T^. 

Cauchy problems for abstract wave equations. Let us consider three Hilbert spaces 
with continuous inclusions 

D{A) CV CH. (84) 

The inner produces and norms of V and H will be recognized with the name of the space 
as a subscript. Let A : D{A) H he a bounded linear operator satisfying: 

(a) an abstract Green's identity 

{Au,v)h + iu,v)v = \/ueD{A) veV (85) 

(b) a surjectivity condition 

I - A: D{A) ^ H is onto, (86) 
/ being the inclusion operator. 
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Then, for arbitrary {uq,Vq) G D{A) x V, the initial value problem 

u = Au t>0, (87a) 
u{0) = Uo, (87b) 
u{0) = vo, (87c) 

has a unique solution 

u e C:2([0, oo); H) n C^([0, oo); V) n C([0, oo); D{A)). (88) 
Moreover, the energy 

e(t) := wmwi + mmi 

is constant as a function of t G [0,oo). Also, the injections f lH^ are dense. The above 
result can be proved by considering the unbounded operator in x if 

: D{A) ^VxH, D{A) := D{A) xV cVxH (89) 

and showing that ±A are maximal dissipative and therefore A is the infinitesimal gener- 
ator of a Cq— group of isometries in the Hilbert space V x H. The possibility of reducing 
all the hypotheses to properties that have to be satisfied by A is related to the fact that 
D{A) is the product space D{A) x V, where the second space is the same as the first 
space in the Hilbert space V x H where the problem is set. 



A 



/ 
A 



Problems displaying rigid motions. Assume that we have three spaces in the same 
conditions above and that in V we also have a seminorm | ■ |y, proceeding from a semi- 
inner product [ ■ , ■ ]v- A space of rigid motions of the system (I57|) is a finite dimensional 
space M such that 

M C D{A), Mcker(v4), 17711^^ = e M, 

and, ii P : H ^ M denotes the orthogonal projection onto M, 

CiWvfy < \v\l + \\Pv\\l < C2\\v\\l WveV. (90) 

In particular, this implies that \v\v = if and only if f G M. The energy of the system 
is now measured in the following form 

e{t):=U^ml + l\u{t)\l. (91) 

If we start with initial conditions (m, 0) G Mx{0}, the solution remains constant over time 
and energy vanishes. If we start with (0,m) G {0} x M, then the solution is u{t) = mt. 
This solution has only kinetic energy and potential energy vanishes identically 

e(t) = i||m||^. 

The conditions are: 
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(modified) abstract Green identity 

{Au,v)h + [u,v]v = WeD{A) veV, (92) 

(b) the surjectivity condition I — A: D{A) H is onto. 

With these hypotheses, problem ( 187|1 has a unique solution with the regularity of (IHHj) for 
any {uq,Vo) G -D(^) x y. Energy, defined with ( 19T|) . is constant in time. 

Remark 8.1. To see how the frame of evolution equations with rigid motions fits into 
the general frame, we need to consider the spaces 

Ho:= {ue H : Pu = 0} = {ue H : {u,m)H = Vm G M}, 

Vo:=VnHo, D{Ao) := D{A) n H^. 

In Vq we consider the norm | ■ |y (see (!90|) ). The operator A has range in Hq, since (p2!) 
implies that 

= (Am, m)/^ + [u, m]y = {Au, m)^ Vm G M. 

Finally, the surjectivity condition implies that I — Aq : D{Aq) H is surjective by the 
same argument. 

Waves in free space. Problem (j4])-([5]), that deals with propagation of initial conditions 
by the wave equation in free space does not fit in the simple frame of this section unless 
we apply some kind of cut-off argument. The main difficulty stems from the definition of 
the energy space. 

(1) It is well known that the closure of the space of smooth compactly supported func- 
tions with the norm ||V ■ Urs is not a subset of L^(M'^). Instead, the resulting 
space can be characterized as a weighted Sobolev space (see [16l [131 [2] ) . The simple 
frame with three spaces (-0(^4) as the domain of the operator, V as the potential 
energy space and H as the space where kinetic energy is measured) cannot be used 
for the corresponding second order Cauchy problem. Instead, the problem has to 
be rewritten as a first order system with four Hilbert spaces involved (compare with 
(!89|) ). Conditions on the operator need to be written in a much more complicated 
form in order to show existence, uniqueness and energy conservation. 

(2) The two dimensional case is even more involved. As explained in [13j, the closure of 
the space of smooth compactly supported functions with || V ■ ||k2 cannot be under- 
stood as a space of functions in any natural way. This adds another complication 
to the four space setting that is needed in the three dimensional case. 

Since the present work deals with compactly supported initial data and speed of propa- 
gation of waves is finite, the strategy of cutting off the analytical domain can be applied 
for any finite time interval, which is enough for our purposes. 
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